Green's Functions and the Adiabatic Hyperspherical Method 
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for a (i-dimensional system. In Section |III| the Green's 
function is applied to the problem of three particles with 



regularized, zero-range, s-wave interactions. Section IV 



applies this result to the three lowest hyperfine states of 
6 Li, and gives a simple description of the scaling behav- 
ior of three-body recombination events that result in trap 
losses. Finally, in Section|V]we summarize the results and 
suggest further avenues of inquiry. 
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THE HYPERANGULAR GREEN'S 
FUNCTION 



Seth T. Rittenhouse, 1 ' 2 N. P. Mehta, 1 - 3 and Chris H. Greene 1 

1 Department of Physics and JILA, University of Colorado, Boulder, CO 80309 

2 ITAMP, Harvard- Smithsonian Center for Astrophysics, Cambridge, MA 02138 

Department of Physics, Grinnell College, Grinnell, I A 50112 

We address the few-body problem using the adiabatic hyperspherical representation. A general 
form for the byperangular Green's function in d-dimensions is derived. The resulting Lippmann- 
Schwinger equation is solved for the case of three-particles with s-wave zero-range interactions. 
Identical particle symmetry is incorporated in a general and intuitive way. Complete semi-analytic 
expressions for the nonadiabatic channel couplings are derived. Finally, a model to describe the 
atom-loss due to three-body recombination for a three-component fermi-gas of 6 Li atoms is pre- 
sented. 



I. INTRODUCTION 

In recent years, there has been extensive theoretical 
and experimental interest in the area of few-body physics, 
most notably in the famous effect first predicted in 1970 
by Vitaly Efimov [1 . Efimov studied a three boson 
system with short-range two-body interactions in which 
each two-body system is infinitesimally close to forming 
a bound state; that is, the s-wave scattering length is 
infinite, or at least very large in magnitude. Quantita- 
tively, Efimov and later others, found that this effect is 
described by a simple wave function in the hyperspher- 
ical representation [5]-[S]. A quantitative understanding 
of three-body scattering [2] [5] has given experiments the 
tools to examine three-body processes in dilute gas sys- 
tems and has led to a wealth of experimental evidence 
for the Efimov effect [THT2"]. More recently, predictions 
relating to the four-body loss rate [T3] have given an- 
other handle on characterizing an Efimov resonance. The 
experimental realization of these predictions swiftly fol- 
lowed 0E|. 

The sticking point of the adiabatic hyperspherical 
method lies in solving the adiabatic Schrodinger equa- 
tion. Often solving this (d — l)-dimensional equation is 
as hard as solving the total d-dimensional Schrodinger 
equation in the first place. Having a variety of methods 
available is therefore helpful. The benefit of using the 
adiabatic hyperspherical method comes from the simple 
final interpretation that can often be applied to the re- 
sulting coupled set of one-dimensional equations in the 
hyperradius 15 . For instance, in the three-body prob- 
lem, if two particles can form a bound state, then one 
of the resulting scattering channels consists of an atom 
and a dimer colliding. In the adiabatic hyperspherical 
method this type of fragmentation channel arises natu- 
rally as one of the discrete solutions to the hyperangu- 
lar equations. In this paper, we derive the hyperangular 
Green's function for an arbitrary d-dimensional system, 
which can then be used in a hyperangular Lippmann- 
Schwinger equation to extract the adiabatic hyperradial 
potential curves. 

This article is arranged as follows, in Section [TT] we de- 
rive a general form of the hyperangular Green's function 



The adiabatic hyperspherical method has proven useful 
for analyzing many few-body systems [TM2"T] . The heart 
of this method lies in treating the overall size of the sys- 



tem, the hyperradius R defined by y/JlR = v/X)j=i Mi 21 ! > 
as an adiabatic parameter. Here pn is the mass scale as- 
sociated with the «th cartesian coordinate and (x is the 
reduced mass associated with the hyperradius. For a sys- 
tem of N particles, Xi denote the d = 3N — 3 cartesian 
components needed to specify the relative positions of the 
N paticles. In this representation, the total wavefunction 
is written: 

* (R, n) = J2 R- (d ' 1)/2 F n (R) $ n (R; n) , (l) 

n 

where the adiabatic eigenfunctions $„(i?;£!) satisfy the 
fixed- R Schrodinger equation: 
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$ n (R;Q,)=u n (R)$ n (R;Q,), (2) 



Here, A is the grand angular momentum operator defined 
by 
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Inserting Eq. (nl into the full (time-independent) 
Schrodinger equation takes a d-dimensional partial dif- 
ferential equation to a set of coupled one-dimensional 



differential equations: 
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F m (i?) = EF n (R) (5) 



The non-adiabatic coupling matrices P and Q in Eq. pi) 
are denned as 



/-',„„ ($ m (R;Cl)\-^-$ n (R;Cl) 
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The integrals in Eqs. ffify and Q are taken only over the 
d — 1 hyperangles collectively denoted 0. 

Approximate solutions can be found by solving the un- 
coupled system of equations, which are referred to as the 
adiabatic approximation: 
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The ground state eigenenergy Eq. ([8]) is a variational up- 
per bound to the exact ground state energy from Eq. ([5]) . 
Another variant of this method is frequently denoted 
the Born-Oppenheimer approximation, with the diagonal 

h 2 

correction Qnn(R) to the potential u n (R) omitted. 

2/j 

These two approximations will be the main focus of this 
paper, while the non-adiabatic couplings P nm (R) will be 
used to describe Landau-Zener-Stueckelberg transitions 
between the different hyperradial channels u n (R). Once 
the adiabatic potentials have been found, much of the 
intuition of simple one dimensional Schrodinger quan- 
tum mechanics can be brought to bear upon the prob- 
lem. Unfortunately, obtaining these potentials can be 
prohibitively difficult in many problems; the development 
of efficient ways to calculate them is desirable. 

This section derives the free space hyperangular 
Green's function for an arbitrary d dimensional space 
such as an A-body system with d = 3 (N — 1) with the 
center of mass coordinate removed. This Green's func- 
tion can then be used to recast equation Eq. pi) into 
an integral Lippmann-Schwinger (LS) equation. The d 
dimensional Laplacian written in hyperspherical coordi- 
nates is given in Ref. [TB] as 

V 2_ 1 d 2 (d _ 1)/2 (d-l)(d-3) A 2 

The hyperangular Green's function is given as the solu- 
tion to 

[A 2 - v (v + d - 2)] G v (SI, CI') = 6 d (Cl~ SI') . (10) 



Here CI stands for the d — 1 hyperangular coordinates 
needed to describe the surface of a d dimensional hy- 
persphere and S d (CI — SI') is the Dirac J-function in the 
hyperangular coordinates, i.e. 5 (CI — SI') = if CI ^ SI' 
and J S d (CI — SI') dSl = 1. The Green's function can be 
found in several forms, including the full hyperspherical 
harmonic expansion [22] . and has been given in closed 
form given by Szmytkowski in Ref. |23j . 

The simplest derivation of the Green's function relies 
on the completeness of hyperspherical harmonics: 

Y,Y^(Sl')Y x ,(Sl)=5 d (Sl-Sl'). (11) 

A/; 

The function Y\^ is the solution to the eigenvalue equa- 
tion 



A 2 y v (Cl) = A (A + d - 2) Y Xfi (SI) 



(12) 



Here A is the hyperangular momentum quantum number, 
and ix enumerates the degenerate states. These functions 
are generally expressed as products of Jacobi polynomi- 
als for any number of dimensions, and are thoroughly 
described by a number of authors (See Refs. [HI [24] for 
some examples) . They are simply an extension of normal 
spherical harmonics to higher dimension. 



Equation ( 11 ) can be used in conjunction with Eq. ( 12 ) 
to find G u (Cl,Sl') P2J: 



G v (Cl,Cl') = Y, 



Y^(Cl')Y^(Sl) 



Xfi 



\(\ + d-2)-v(v + d-2)' 



(13) 



Unfortunately, eigenfunction expansions of Green's func- 
tion often have slow convergence with respect to the num- 
ber of eigenfunctions, making them unsuitable for numer- 
ical calculations. The closed form of the Green's function 
from Ref. [23] is given as 



G v (Cl, SI') 



(d — 2) Sd sin7w 



Cj.dr.2y2 (^_ R . ft\ ; (u) 



where C" is a Gegenbauer function, Sd is the surface 
area of the d-dimensional unit hypersphere: Sd — J dSl — 
27r d / 2 /r (d/2), and R • R! is the cosine of the angle be- 
tween the two normalized hypervectors R and R' . Here 
v is defined by Eq. (10). While Eq. 14 has a pleasing, 
compact form it is often divergent at critical points. For 
instance if v is non-integer valued, then G v (Cl, Cl') di- 
verges as R ■ R! —* 1 . 

For these reasons, it is convenient to find a third form 
of the Green's function. The first step in this derivation 
relies on the division of the total d dimensional space into 
two subspaces. For the purposes of this work, we will 
assume that the dimension of the two subspaces are both 
greater than 2: i.e. di,d^ > 2. The two subspaces are 
each described by sub-hyperspherical coordinates. The 
the two resulting sub-hyperradii can then be related to 



the total hyperradius as 



R± = i?sina, 
i? 2 = -Rcosa, 

< a < it/2. 



(15) 



Reference |24| details how the hyperangular momentum 
can be written in terms of the sub- hyperangular momenta 
as in Eq. (17 1. With the following definitions: 



A?y Al/11 (Oi) = Xi (Ai + di - 2) y AlM1 (fii) 
A^r A2A12 (o 2 ) = a 2 (A 2 + d 2 - 2) r A2M2 (fi 2 ) 



(16) 



the Green's function can be expanded using the com- 
pleteness of the sub-hyperspherical harmonics [viz. 
Eq. (18)]. Substituting the expansion Eq. (18) into 
we find that the latter is satisfied if and only if 
is satisfied; S (a — a') is a Dirac <5-function and 



Eq 



Eq. (19 



10 



the denominator on the left hand side of Eq. ( 19 1 arises 



from the hyperangular volume element associated with 
the angle a (See Refs. [TUH2] f° r details). 
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The general one-dimensional Green's function for any differential equation of Sturm-Liouville form Eq. ( 19 ) is: 

-/aWOO/a^O^) 



a 1. do. { !\ 



{sinaf 1 ' 1 (cos a) 1 * 2 ' 1 W [/+, /" 



(20) 



where W[f + ,f ] = f + f ' — f f +l is the Wronskian [25] and a<(>) = min(a, a 1 ) (max (a, a')). The functions 
/ + (a) and /~ (a) are regular at a = and a = 7r/2 respectively and satisfy the homogeneous version of Eq. (19 1. 
The solutions / + (a) and /~ (a) are given in Ref. [5S] as 



/S 2I/ (a) = (sin Al acos A2 a) 2 F 1 ( 



A1+A2-1; v+\ 1 +\ 2 +d-2 . 2\±+d± _ l=pcos2a 



UAiA 2 ^' * AiA 2 vJ / . 



2p /' 2Ai+di \ p /- 2A 2 +d 2 \ 
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(21) 
(22) 



where 2 -F\ (a, 6;c, x) is a hypergeometric function, A + 
Ai, d+ = di, A_ = A 2 , and d_ = d 2 . 



III. THE THREE-BODY PROBLEM WITH 
ZERO-RANGE INTERACTIONS 



the three body problem with regularized, zero-range, s- 
wave, pseudo-potential interactions. This problem has 
been well studied by a variety of sources 4, 27 -29J. The 
full Hamiltonian for the untrapped system is given by 



In this section we show the utility of the Green's func- 
tion developed in the previous section by applying it to 



^ = E-2^7 V ' + E%(^)> ( 23 ) 

i=l - ' i>j 



where f*j is the position of the ith particle, and V 2 is the 
Laplacian for fi . The interaction is given by 



Vij (r«) 



47r/l 2 , 

2//j 
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l3J dri/ v ' 



(24) 



where ay is the s-wave scattering length between par- 
ticles i and j and (iij is the two body reduced mass, 
fiij = miTrij / (mi + rrij) . The pseudo-potential defined 
in this way applies the Bethe-Peierls boundary condi- 
tion to the two-body wave function as r — > 0, i.e. 
-0 (r) — » (1 — dij/r) C for some constant C 30. The 
center of mass can be removed from this system by con- 
verting to a system of Jacobi vectors. Jacobi vectors are 
created for this system by considering the separation vec- 
tor between two of the three particles and then a second 
vector from the center of mass of that two body system to 
the third. The final vector is then just the center of mass 
coordinate. The choice of Jacobi vectors is not unique. 
Here we will need to consider three different Jacobi coor- 
dinate parametrizations each of which is convenient for 
describing one of the three possible two-body interactions 
V (rij ) . In the "odd-man-out" notation these are given 
by 

Pi ] = (n-r^/dk, 



Jk) , I nun + mjTj 

A = d k [ ; — — - r k 



171,1 + TTlj 
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2 _ (m k /n){mi +rrij) 
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where fi is the three-body reduced mass: 
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mim 2 m 3 
nil + m 2 + m 3 



(26) 



The total Hamiltonian can be rewritten in terms of the 
Jacobi coordinates and the center of mass as 



H, 



Htot — H + Hcm, 

-h 2 

2M 



(27) 
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Transforming the Jacobi coordinate piece of the Hamil- 
tonian in Eq. (27) into hyperspherical coordinates using 
Eqs. ^ and |25j yields 



F = --4^ii? 5/2 l5h 
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To apply the adiabatic hyperspherical formulation, the 
hyperangular adiabatic Schrodinger equation must be 



solved: 



Aa + ^£ Vtf ( dfcp «)_, (( , + 4) 



$(i?;fi) =0. 

i<j 

(29) 
This can now be accomplished with the use of the hy- 



perangular Green's function, Eq. (18), in the Lippmann- 
Schwingcr (LS) equation, 



$ (R; Q) = -^p- / dVl'G" (fi, n') 
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(30) 
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where p{ is the fcth Jacobi vector parametrized by 
{R, fi'}. Because the system has been constrained to 
have a constant hyperradius, this is effectively a bound 
state problem; note that Eq. p0| ) has been assumed here 
to have no noninteracting solution at the chosen value 
of v . The hyperradial Hamiltonian from Eq. ([5]) in the 
absence of the non-adiabatic couplings P and Q is given 
in terms of the hyperangular eigenvalue v as 



H B = 



-h 2 d 2 
2/7 aR 2 



U n (R) , 



(31) 



U„ (R) = 
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R 2 



^Cnn \-*v 



To evaluate the integrals in the LS equation, the 
Green's function from Eq. ( 18 ) is expressed in terms of 



the appropriate Jacobi coordinate set for each interaction 
term in the sum, with the hyperangles defined as 



Q« = {u[ 
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where uj\ represent the spherical polar angular coordi- 
nates for (h_- The remaining hyperangle aS ' is defined 



as in Eq. (15), i.e 



pf } =i?sina (fe) , 



(33) 



!h 



i k) =Rcosa^. 



With this choice of hyperangles, it is clear that d\ = d 2 = 
3 and the hyperspherical sub- harmonics Y} (ftj) in 



Eq. ( 18 ) reduce to normal spherical harmonics VLtMi (wj). 
The (5-function implies that the Bethe-Peierls two-body 
boundary condition for each two-body interaction can be 
applied and the third particle can be considered to be far 
away, i.e. 



pi» 



lim $ (R; fl) = 
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Here x/lm is a spherical harmonic describing the free 
space behavior in uj 2 and it carries the total angular 
momentum of the system. The superscript fc again indi- 
cates the odd man out notation. This gives the values of 
the sub-hyperangular momentum quantum numbers in 



the fc Jacobi coordinate system as Ai = and A2 = L, 
which accounts for the s-wave interaction and the total 



angular momentum L. Inserting Eq. (34 1 into Eq. (30) 
gives the hyperangular eigenf unction, 
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where (j,k is the two-body reduced mass labeled in the odd man out notation and the orthonormality of spherical 

(fcV (k)/ . c- J I 

harmonics has been used to evaluate the u\ and UI2 integrals. The d-function in Eq. (24) implies that the integral 



in a^ k '' can be accomplished by evaluating at a 



(fc) 



*( fe )' = 0. 



The analytic equation for the hyperangular eigenfunction in Eq. ( 35 1 is not very useful without knowing the 



hyperangular eigenvalue v(R). To obtain an equation for v the boundary condition given in Eq. (34) must be applied 
again, i.e. 
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To evaluate the limit on the right hand side of this, we 

must determine the values of the k 7^ k' Jacobi coordi- 

(V) J 

nates in the limit p\ — > 0. Equations (25) and (33) 

give, for k 7^ fc', 



lim a<*> = kk , -- 

lim 4 « ~A 
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(37) 
(38) 



The hyperangular eigenvalue, v, is found by solving the 
closed form transcendental equation, 



det (M - 



0, 



(41) 



for any given total angular momentum L, any set of s- 
wave scattering lengths a\ k > and arbitrary masses. 



A. Imposing symmetry 



Note that if / is regular at fikk' , then 
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Using this and evaluating the limits in Eq. ( 36 1 yields a 
matrix equation for C L lj\ 
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The hyperangular eigenvalues for the general three- 
body problem with arbitrary exchange symmetry can 
be found by solving the transcendental equation (41), 



but the system can be simplified by considering differ- 
ent permutation symmetries and imposing those symme- 

(fc) 
tries on the boundary conditions C^ M . For example, if 

the particles in question are identical bosons, permuta- 
tion cannot have any effect on the wave function. Thus, 
if two particles are exchanged in the two-body subsys- 
tem, the boundary condition must remain the same, i.e. 
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A complete list of the possible exchange symmetries is 
given in Table [TJ 

To illustrate this post-symmetrization, we apply the 
identical boson symmetry with L = to Eq. ( 40 ) re- 



suiting in the well known transcendental equation for v 
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In the limit where R -C \a\ the first solution to this tran- 
scendental equation gives a super-critical attractive 1/R 2 
effective potential, 



U(R) 



1/4 



2[i R 2 
s = 1.00624. 



(43) 



This attractive potential is the source of the famous 
Efimov effect, where an effective attractive dipole-type 
potential supports an infinite set of three-body bound 
states that accumulate at the non-interacting three-body 
threshold, E = 0. 
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TABLE I: The possible permutation symmetries that may be 
imposed on the three body system with s-wave interactions 
are given with the appropriate boundary conditions. B stands 
for a boson, F for a fermion and X for a distinguishable par- 
ticle with an arbitrary mass. 



B. Non-adiabatic couplings 

As in any adiabatic treatment, the effective hyperra- 
dial potentials are coupled by non-adiabatic terms that 
arise from the hyperradial dependence of the hyperan- 
gular channel functions. These couplings come in the 
form of the P- and Q-matrices in Eq. ^\. To find the 
non-adiabatic coupling matrices, we apply the methods 
of Ref. [33] . The details of the derivation are shown in 
Appendix A, the result of which gives the semi-analytic 
expressions for the matrix elements P mn '- 
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for notational simplicity, 



we have set e„ — 



Here, 

(v n + 2) z and all primes indicate a derivative with re- 
spect to R (e.g. e' n = de n /dR). Because the hyperan- 
gular eigenfunctions are orthonormal, the diagonal part 



of the P matrix is zero, i.e. P„„ = 



ld_ 
2dR 



($„!$„) = 0. 



Equation (45) gives the normalization condition for $ n , 



with an overall phase that is free. This overall phase is 
chosen here so that ^ fc C„ is positive. A similar deriva- 
tion provides the matrix elements Q m n'- 
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When the symmetries given in Table I are used, there 
can be a considerable simplification of the expressions for 
P mn and Q mn - For a system of identical bosons where 

fl (l) = ft (2) = (3) = a , dl = d2 = d3 
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,1, Pmn and Qmn are given by 
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which are in agreement with previously calculated nona- 
diabatic corrections for the three identical boson system 

123. 



IV. THREE DISTINGUISHABLE 
INTERACTING PARTICLES. 



In this section the adiabatic three-body potentials and 
the non-adiabatic couplings are applied to the case of 
three distinguishable equal-mass particles. This system 
has been realized, for instance, in ultracold three compo- 
nent Fermi gases of 6 Li atoms [H [10] which has sparked 
a great deal of recent theoretical interest [33-36 . The 
scaling behaviors and recombination rates we discuss in 
this section can be found in Ref. [35]. We derive them 
here to illustrate the power of the methods presented in 
this paper. The scattering lengths near the resonance 
positions used here, as functions of magnetic field, are 
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FIG. 1: (color online) All possible s-wave scattering lengths 
are shown for the lowest 3 Zeeman states of Li 6 from Ref. 
|37| . Each marked region gives a different set of length scale 
discrepancies. Here cr*' is the scattering length between two 
atoms in states |i) and \j) with k as the component not in- 
volved in the interaction. 



? iven in Refs. [101137] by 

A 



a« - a b 



1 



[1 + a (B - B Q )} , (48) 



B-B a 

for k = 1 : 

a b = -1450a , B Q = 834.15 G, 

A = 300 G and a = 4 x 10~ 4 G _1 ; 

for k = 2 : 

a b = -1727a , B = 690.4 G 

A = 122.2 G and a = 2 x 10~ 4 G _1 ; 

for k = 3 : 

a b = -1490a , B = 811.22 G, 

A = 222.3 G and a = 3.95 x 10~ 4 G _1 ; 

where clq is the Bohr radius. The Fano-Feshbach reso- 
nances in this system allow for a large variety of tunable 
interactions. 

This series of overlapping resonances produces five dif- 
ferent regions of magnetic field, shown in Fig. [I] near 
the three resonance positions, each possessing distinct 
behavior. In all five regions the scattering lengths are 
much larger than the effective range, allowing for the use 
of the zero-range interaction assumptions. Table [Tl] shows 
the various length scale disparities in these regions. 

Figure [2Fa) shows an example of the lowest four hy- 
perangular eigenvalues (y + 2) obtained from solving 
Eq. (J40]) for aW = a ^ and a^ = 1000aW. This is pro- 
vided as an example that is qualitatively similar to the 
behavior of the system in region I. When the hyperra- 
dius is in a region where all other length scales are much 
different, the hyperangular eigenvalue (y + 2) becomes 
constant, or, in the case of 2-body bound states, becomes 



Region 




I 


ro « a (3) < a (1) « a< 2 > 


a^,a^,a^ >0 


II 


r < a (3) ~ a (1) < 


a< 2 > 




a^ <0;a^,a^>0 


III 


r < 


a< 2 > 


«a< 1 »,a( 3 ) 


a^ <0;a^,a^>0 


IV 


r < 


a< 2 > 


< 


a« 


,a (3) 


a( 2 >,a«<0;a( 3 >>0 


V 


r < 


a< 2 > 


< 


a W 


< 


a< 3 > 


aW,a( 2 », a < 3 » <0 



TABLE II: The possible tunable interaction regimes near the 
resonances of 6 Li are given. 
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FIG. 2: (color online) (a) For an example system having 
(V- ' = a^ 3 ' and a' 2 ' = 1000a' 1 , the first four hyperangular 
eigenvalues are shown as functions of the hyperradius. The 
solid black horizontal lines show the expected behavior for 
3 identical resonantly interacting bosons. The dashed line 
gives the behavior of two identical fermions interacting reso- 
nantly with a third distinguishable particle. Dotted lines give 
the expected universal behavior for a single resonant scat- 
tering length. Finally, the dot-dashed line is the lowest ex- 
pected free space eigenvalue for three distinguishable free par- 
ticles. (b)The coupling strengths between the third and fourth 
(purple solid curve), the first and fourth (red dashed curve), 
and the first and third (black dotted curve) adiabatic poten- 
tials arc shown as a function of R. 



proportional to R 2 . This behavior can be interpreted as 
giving a universal set of potential curves from Eq. (31 1 



For example in region I where r$ <C a 1 - 3 ' < a*- 1 -* <C a [ ' 
there are three hyperradial regions: r <C R <C a^ < 
a (i) < a (2). r < a (3) < a (i) < i? < a (2). an d 



ro 
(3) < (i) 



r ^C a w < a (1> <C a^ <C R- In each region the hy- 
perangular eigenvalues take on the universal value that 
is expected for resonant interactions [H (57J [55] . 

Figure [3] schematically shows the behavior of the first 
few hyperradial effective potentials from Eq. (31). The 



grey areas are the regions where potentials are transition- 
ing from one universal behavior to the next. The zero- 
range pseudo-potential cannot describe the short range 
details of the interaction, meaning that the potentials 
found here are only valid for fi> tq where rrj is a short 



range parameter shown schematically as the labeled blue 
region of Fig. [3j 

r a^\a^ a ^ 



ft? 




FIG. 3: (color online) A schematic picture of the first four 
hyperradial potentials in region I is shown. The grey areas, 
labeled by the appropriate scattering lengths, indicate regions 
where the potentials are changing from one universal behavior 
to another. The blue region, labeled by ro, indicates the short 
range region where the zero-range pseudo-potential not longer 
can be applied. 

Figure [2Vb) shows the coupling strength, 
P^ n /2fj,[u m \R) — u n (12)], between the different 
potentials. The places where this coupling peaks are 
the points where a transition between curves is the 
most probable. Figures |i[a-e) are examples of the 
hyperangular eigenvalues found in each region. The 
magnetic field at which each set of eigenvalues are found 
is shown as dotted lines in Fig. BFf) from left to right for 
Fig. WTa-e) respectively. In each figure the hyperangular 
eigenvalue can be seen flattening out to a universal 
constant in each region of length scale discrepancy. As 
the magnetic field is scanned through each resonance, 
one two-body bound state becomes a virtual state. This 
behavior can be seen in the hyperangular eigenvalues 
that diverge toward — oo. As each resonance is crossed, 
one of the hyperangular eigenvalue curves goes from 
diverging to — oo to converging to (v + 2) — > 4. 

As a final examination of this system, we extract the 
scaling of the low energy three body recombination rate, 
i.e. the rate at which three particles collide and form 
a dimer and a free particle. The lowest 3-body curve, 
the lowest potential that goes to the three-free-particle 
threshold, is the potential that dominates this process. 
Contributions from higher hyperradial potentials will be 
suppressed due to larger tunneling barriers. One limi- 
tation of the zero-range pseudo-potential is that it only 
admits at most one dimer of each type. The process of 
three-body recombination releases the binding energy of 
the dimer state as kinetic energy between the dimer and 



FIG. 4: (color online) (a)-(e) Examples of the hyperangu- 
lar eigenvalues from each region of magnetic field are shown 
as a function of the hyperradius in atomic units, (f) The 
three s-wave scattering lengths are shown as a reference plot- 
ted versus the magnetic field strength. The vertical dot- 
ted lines, from right to left, show the magnetic fields at 
which the hyperangular eigenvalues in (a)-(e) were found, 
B = 675, 695, 805, 820 and 845 gauss respectively. 



remaining particle. For the purposes of this study we will 
concentrate on the three-body recombination processes 
that result in trap loss processes, where the energy re- 
leased in the recombination can be assumed sufficient to 
eject the remaining fragments from a trap. 

The event rate coefficient for TV initially unbound par- 
ticles with total orbital angular momentum L to make 
a transition from a hyperspherical potential curve with 
hyperangular eigenvalue A to a lower lying final state is 
given bv [41138] 



K N = 



hk 



-N s 



d-l 



27T rf / 2 



i,f 



(2L + 1) \T fi 



where d is the total dimension of the system (in the case 
of three-body recombination d = 6), T/j = Sfi — Sfi is 
the transition matrix element between an initial three- 
body entrance channel i and a final exit channel /, and 
k = y/2pE/H is the wave number of the asymptotic hy- 
perradial wavefunction. The sum in this equation runs 
over all the initial, asymptotic channels with total angu- 
lar momentum L that contribute to the scattering pro- 
cess. In Eq. (49 1, Ns is the number of permutational 



symmetries in the system. For three distinguishable par- 
ticles Ns — 1, but it can be different, for instance for A" 
identical bosons, Ns = N\. 

In the low energy regime, only the lowest L = initial 
three-body channel will contribute, while higher chan- 
nels will be suppressed. The sum over final T matrix ele- 
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merits can be approximated using the Wentzel-Kramers- 
Brillouin (WKB) phase in the entrance channel [28l [38] : 



£|7>;l 

/ 



= -27 



sinh 2r/ 



2 cos 2 (j> + sinh r\ 



(50) 



where r\ is an imaginary phase which parameterizes the 
losses from the incoming channel. In Eq. (50) 7 is the 



total WKB tunneling integral between the outer classical 
turning point and the hyperradial position at which the 
transition to the outgoing state occurs, i.e. 



7 = Re 



2/x 

F 



[U(R)-E] 



1 
4R 2 " 



dR 



(51) 



where E is the initial three body energy, Rj< is the outer 
classical turning point and Rq is the position at which 
the coupling between the incoming and outgoing channels 
peaks. In Eq. (50) is the WKB phase acumulated in 



any inner attractive well: 



Im 



Rt 



ri„ 



2/t 



[U{R)-E] 



4i? 2 



dR 



(52) 



The extra repulsive 1/4R 2 term in Eqs. (51) and (52) is 
due to the Langer correction [55] . The total T- matrix 
element will depend on the detailed nature of the real 
short range interactions and the behavior of the outgo- 
ing channels, but the scaling behavior with the scattering 
lengths will be determined by Eq. (50). In each region of 



magnetic field, there are different length scale discrepan- 
cies and different numbers of bound states. As a result, 
we will examine each region separately. 



Region I (a. 



a< 3 > « a< 2 >) 



FigureEla) shows the behavior of the first few hyperan- 
gular eigenvalues in region I. The first three eigenvalues 
correspond to dimer states, while the fourth corresponds 
to the lowest three-body potential and is the entrance 
channel that will control three-body recombination. The 
lower two dimer states are relatively deeply bound with 
binding energies, h 2 /ma 2 , on the order of 10 -12 Hartree. 
This is comparable to the trap depth energy of a normal 
magneto-optical trap for experiments with 6 Li HO] , 
meaning that recombination into these dimer channels 
typically releases enough energy to eject the remaining 
dimer-atom system from the trap. 

In the limit where R ^$> a^ 3 \ the three atoms are far 
enough apart to be in the non-interacting regime. This 
means that the hyperangular eigenfunction limits to the 
lowest allowed three-body hyperspherical harmonic with 
its corresponding eigenvalue, (y + 2) — ¥ 4. In this limit 
the hyperradial potential becomes 



U (R-»a 



(3) 



h 2 4 - 1/4 
27 R 2 



(53) 



For very low energy scattering, the classical turning point 
in Eq. ( 50 1 , is approximately 



Rt = 



1 
2k' 



(54) 



In fact, this will be the turning point for all of the three- 
body recombination processes discussed in this section. 

It is possible for recombination to occur directly be- 
tween the lowest three-body curve and the deep dimer 
channels, but this direct process is strongly suppressed 
due to the large tunneling barrier in the three-body po- 
tential at small R. The favored path is through a transi- 
tion to the weakly bound dimer channel, shown schemat- 
ically in Fig. [5] The coupling between the lowest three- 

r aV,aW (a) 







FIG. 5: (Color online) A schematic of the path for three-body 
recombination in region I is shown. Transition regions are 
labeled by the appropriate length scale, and the short-range 
non-universal region is labeled by tq. 

body channel and the weakly bound dimer channel peaks 
at approximately R ~ a} 2 \ while the coupling peak be- 
tween the weakly bound dimer channel and the remain- 
ing two dimer channels occurs at approximately R ~ a*- 3 ' 
~ a^\ In the regime where a^ ~ a^ 3 ' <C R <C af^> the 
three particles are so far apart that they cannot see the 
smaller scattering lengths aP^ and a^ 3 \ but the third 
scattering length is so large compared to the hyperradius 
that it might as well be infinite. This leads to a universal 
potential wh ose hyperangular eigenvalue can be found by 
solving Eq. (40) with a^ = a^ = and a^ -> 00, i.e 



U(R) = 



h 2 p\ 



1/4 



Pi 



2/y 

iy- 



R 2 

2) = 1. 



(55) 



This intermediate universal behavior can clearly be seen 
in Fig. (|2ji). 

The behavior of each channel can be approximated 
by the universal behavior of the hyperradial potential 
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in each region. Under this assumption, using Eq. (50 1, 
the tunneling probability is given by 



-2t 

e ' oc exp 



p'i 



iii 



- 2/ 1(3) dR W~f E 



(56) 



2 I dR 

1.( 2 > 



r 2 h 2 



If the scattering energy is very small, E <C h 2 /m [a^] , 
then the energy dependence in these integrands becomes 
negligible leaving, 



-27 



(a^a^j 



(57) 



Inserting this in for the T- matrix element in Eq. ( 49 ) 



gives the scaling behavior of the recombination rate with 
the scattering lengths [28 : 



^ 3 °c(>v 3 >y 



(58) 



It was assumed here the final transition occurs at R ~ 
a^ leading to the scaling behavior with a^ 3 \ but the 
transition could just as easily have occurred at R ~ oA 1 '. 
a^ 1 ' and a^ 3 ' are approximately equal here, and which 
one dominates the transition depends on the short range 
behavior of the real two-body interaction. To extract 
the scaling behavior with res pec t to a' 1 - 1 , one can simply 
replace a' 3 -* with a^ in Eq. (58) as long as a^ 1 ' and a^ 
are approximately equal. 



2. Region II (a 



(i) 



a (3) « 



,(2) 



The recombination in region II is simpler than in re- 
gion I as there is no weakly bound intermediate state. 
Again, we assume that the trap loss recombination is 
dominated by transitions to the two remaining dimer 
states seen in Fig. Kb) . The lowest three-body potential 
has coupling to these channels that peaks at R ~ a' 1 - 1 
and R ~ a^ 3 \ For R 3> \a^\ the hyperangular eigen- 
value takes on the non-interacting value [y + 2) — > 2. 
For a^\a^ <C R <C M 2 ''! the universal hyperangular 
eigenvalue (v + 2) = p\ = 1 is seen again [2"TH2"5] , Ignor- 
ing the transitional region between these two regimes the 
transition probability is given by 



e 2l oc exp 




I, the scaling behavior with respect to a^ can be found 



by simply replacing a^ with a' 1 - 1 in Eq. (58) as long as 
a^ 1 ) and a^ are close. 



3. Region III ( a (2) < a (1) ~ a (3) ) and Region IV 

(]«( 2 »|«|fl< 1 )]~fll 3 >) 

In Region III, none of the dimers predicted by the zero- 
range model have enough binding energy to cause trap 
loss. While recombination can occur into these channels, 
we will focus on the process of recombination to deeply 
bound states here. In reality, the deep interaction poten- 
tial between two Li atom in different spin states admits 
many deeply bound dimer states, and a true hyperspher- 
ical description of the system would have channels going 
to each possible dimer-atom threshold. The energy re- 
leased in recombining into these deep states is enough 
to kick the atoms out of any normal trap. Because the 
deeply bound states are of the size of the range of the in- 
teraction, coupling to the deeply bound hyperradial chan- 
nels will peak at small hyperradius, R ~ r , and the rate 
can be found by studying the tunneling probability of 
reaching these states. 

As with the recombination process in Region I, the 
most favorable pathway involves multiple steps. Starting 
from the lowest three-body channel, a transition is made 
to either the first or second weakly bound dimer chan- 
nel. Because a*- 1 ' and a^ 3 ' are similar in magnitude, the 
coupling to these channels peaks in the same region. If 
the transition is made to the highest dimer channel, then 
another transition is made directly to the second. 

This pathway is shown schematically in Fig. [6] An 
interesting thing occurs in the lowest weakly bound po- 
tential when la' 2 ' I <C R <§C a^ ~ a^: the universal 
potential becomes attractive. This region of attractive 
potential gives rise to a number of phenomena. For in- 
stance, in the limit a^ x ',a^ 3 ' — > oo, the universal attrac- 
tive potential supports an infinite number of geometri- 
cally spaced three-body bound states, giving rise to the 
Efimov effect. In the process of three-body recombina- 
tion to deeply-bound dimer states, though, there is no 
tunnelling suppression in this channel, and the hyperra- 
dial wavefunction merely accumulates phase in this re- 
gion. As a result the WKB tunnelling probability is con- 



trolled by the transition at R 



,,(!). /i(3). 



e 2j oc exp 



R T 


1 4 


dR\ 

(i) V 


'r 2 



(61) 



Inserting this into Eq. p) gives a recombination rate Inserting this into Eq. (|49j) gives the scaling of three-body 



that has the same scaling behavior as in region I 



recombination to deep dimer states as 



Ko a 



(a( 2 V 3 >)' 



(60) 



i^oc 



,.(i) 



(62) 



Again it is assumed that the final transition occurs at 
R ~ aS 3 \ but it could occur at a^ 1 ' as well. As in Region 



Again, it is assumed here that aP^ and a' 3 -* are similar in 
magnitude. If this is not the case, for instance if a^ 3> 
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FIG. 6: (color online) A schematic of the potentials and the 
path for three-body recombination in region III is shown. 
Again the labeled grey regions indicate a transition from one 
universal potential behavior to another. The green line rep- 
resents the hyperradial potential for a deeply bound dimer 
state. The blue area labeled by ro is the short range region 
not described by zero-range interactions. 



a^ , then a scaling behavior similar to that of Eq. ( 58 ) 
is recovered: 



Kocx 



aWa® 



(63) 



In Region IV there is only a single weekly bound dimer 
state available, and trap loss will occur through recombi- 
nation to deeply bound dimers. The path here is similar 
to that of Region III, where a transition happens from 
the lowest three-body channel to the weakly bound dimer 
channel. From there the hyperradial wavefunction can go 
to the small R region without further suppression. This 
process then yields the sam e three-body recombination 
scaling behavior as Eq. (62) when a^ ~ | a ^|- When 
a^ 3> \a^\ , the scaling predicted by Eq. (63) is recov- 
ered. 



4. Region V ( \a ( >\ < 



,(3) 



«aW) 



In this regime the recombination process is entirely 
controlled by the lowest three-body channel, shown 
schematically in Fig. [7] The hyperradial potential has 
three universal regimes. The first, when r <C R <C 
|a (2) | < |a (3) | < I *" 1 '!; is identical to that of three 
strongly interacting bosons. The hyperangular eigen- 



value, (is + 2) , is the first solution to Eq. (42) in the 
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FIG. 7: (Color online) A schematic of the lowest hyperradial 
potential is shown with the path for three-body recombination 
to deeply bound states. The green line represents the hyper- 
radial potential for a deeply bound dimer state. Labeled grey 
areas indicate transition regions from one universal behav- 
ior to another, and the blue region indicates the short range 
regime. 



limit where R/a — > 0, yielding the hyperradial potential, 



U(R) = 



h 2 - ( SQ ) 2 - 1/4 



1[i R 2 

s = 1.0062. 

In the next regime, when tq -C 



(64) 



,(2) I 



<C R < 



,(3) I 



< 



>(1) 



, the three particles are far enough apart so as not 
to see the smallest scattering length. As a result the 



hyperangular eigenvalue is governed by Eq. ( 40 1 with the 
BBX symmetry of Table |T] imposed: 



U(R) = 



( Sl ) 2 -l/4 



2[i R 2 

si = 0.4137. 



(65) 



In the regime where ro <C \a^\ <§C \a^\ Cfi<C I - 1 ' 1 '! 
there is only one scattering length that is seen by the sys- 
tem, and the universal potential becomes that of Eq. ( 55 1 



In the final regime, where the hyperradius is much larger 
than all of the scattering lengths, the potential goes to the 
non-interacting behavior of a hyperspherical harmonic. 

The transition to a deeply bound dimer state occurs 
at R ~ ?' following the path shown in Fig. [7J To get to 
this region, the wavefunction must first tunnel through a 
barrier, leading to suppression of the recombination rate. 
Once through the barrier, the wavefunction accumulates 
phase in the attractive potential regime. If enough phase 
can be accumulated in this regime, then a three-body 
bound state (a so called Efimov state) can be present 
leading to a resonance in the recombination rate. The 
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final recombination rate for this process is [H [S7J I2H] 

sinh 2rj 



K 3 (xA 



cos 2 {4>wkb) + sinh 77 



(66) 



where r\ is controlled by the short range properties of the 
system, 4>wkb is the WKB phase accumulated in the 
attractive regime from r$ to I a^ I , and A is proportional 
to the tunneling suppression through the barrier: 



Aoc 
^wkb = si In 



a^a^ 



a( 2 ) 



s In 




(67) 
(68) 



Figure [8] (a) show the log of the three-body recombina- 
tion rate in arbitrary units as a function of the magnetic 
field with 77 = 0.001. The short range length scale here is 
chosen to be approximately the van der Waals length of 
6 Li, tq — rd w 30 atomic units. Figure M (b) shows the 
scattering lengths in the same region of magnetic fields 
for reference. An Efimov resonance can clearly be seen at 
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FIG. 8: (c olor online) (a)The three-body recombination rate 
from Eq. (661 for 6 Li is shown in arbitrary units as a func- 
tion of magnetic field with 77 = 0.001 and the short range 
length scale chosen to be approximately the van der Waals 
length, ro = rd ~ 30 a.u. The large y-axis tick marks indi- 
cate orders of maginitude. (b) The three scattering lengths 
a' 1 ' (solid black curve), a> 2 ' (dashed red curve) and a' 3 ' (dot- 
dashed blue curve) are shown in atomic units as a function of 
magnetic field in Region V. 

B = 942 G when 4>wkb = 37r/2. This is in rough agree- 
ment with the predicted position of B = 1160 G found in 
Ref. [33J . The exact position of this resonance is some- 
what sensitive to the short-range length scale r$ which 
should be fit to experimental data. We have chosen ro as 
the Van der Waals length here for illustrative purposes. 
A WKB phase is 3tt/2 at the resonance indicates that 



this corresponds to the second Efimov state intersecting 
the continuum. The first Efimov state remains bound 
throughout this region. Because a^ becomes resonantly 

large asB-> 834.15 G, the [ a ( 3 )] 2 scaling from Eq. (67) 
gives the large recombination rate seen in the lower field 
region of Fig. M (a). 

With three overlapping resonances, 6 Li provides a rich 
hunting ground for the study of three-body physics. Fur- 
ther, because it is a Fermionic atom, three-body interac- 
tions involving only two of the three lowest components 
are strongly suppressed meaning that the majority of the 
three-body physics is controlled by a system of three dis- 
tinguishable particles. While only the processes of three- 
body recombination that lead to trap losses were studied 
in this section, there is still a rich and complex array of 
behaviors not discussed that can be described using the 
model presented here. 



V. SUMMARY 

In this work we have developed a new form for the 
hyperangular Green's function in arbitrary dimensions. 
The derivation of the Green's function is simple and fol- 
lows easily from a standard Sturm-Liouville problem. By 
dividing a d dimensional space into physically meaning- 
ful subspaces, this new Green's function avoids the slow 
convergence often seen in spectral expansions form, while 
maintaining a physically intuitive set of hyperangular co- 
ordinates. 

We have also used the hyperangular Green's func- 
tion to solve the three-body problem with zero-range s- 
wave interaction for arbitrary scattering lengths, parti- 
cle masses and total angular momentum. With simple 
root finding, the adiabatic hyperangular channel func- 
tions and adiabatic potentials can be extracted. The 
resulting transcendental equation is in exact agreement 
with that derived using Fadeev like decompositions. To 
complete the problem, we have also derived, for the first 
time, general expressions for the non-adiabatic correc- 
tions to the potentials that are analytic up to root find- 
ing. 

The results of the general three-body problem were 
then applied to the three lowest hyperfine components 
of 6 Li near a set of overlapping resonances. By a simple 
WKB formalism, the scaling behavior of rate constant for 
trap loss three-body recombination events was extracted 
throughout the overlapping two-body resonances. Sig- 
natures of an Efimov style resonance are also predicted 
to appear at high field strengths. Throughout the reso- 
nances, all of the scattering lengths are very large com- 
pared to the length scale of the two-body interaction, 
indicating that the results presented here are univer- 
sal. The simple and intuitive nature of the Lippmann- 
Schwinger equation in the three-body problem indicates 
that this Green's function based method may be appli- 
cable in the context of the four-body problem, but this 
extension is the subject of ongoing inquiry. 
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Appendix A 

In this appendix we sketch the derivation of the for- 
mulas for the non-adiabatic P and Q matrix elements 
given in Eqs. (44) and (46). We begin by considering 



matrix elements dealing with the derivative of the Adia- 
batic Schrodinger equation: 



-4($ m |$„; 



«|(A 2 -e m 
-($ m |(A 2 -e, ; 



I K) = o 



(69) 



where e n — v n (v n + 4) is the hyperangular eigenvalue of 
the nth adiabatic eigenf unction, and the prime indicates 
a hyperradial derivative has been taken. Taking the dif- 
ference of these leads to an equation for the non-adiabatic 
coupling matrix element P mn for m =/= n: 

« |A 2 | $ m )-($ m |A 2 | $'„)-(£,„ - £„) Pmn+SmnS'n = 0- 

(70) 
The difference ($'„ |A 2 | $ m ) - ($ m |A 2 | $' n ) is given by 
the boundary conditions of the wave functions <!>,„ and 
$„ at the coalescence points: 
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«|A 2 |$ ro >-($ m |A 2 |$;> = ^ 



a (k) r JLrW - r( k )JL f a(k) n(k) 
dtR m dR n " ! ■ " " 



9i? V d k R 






(fc) 



4-R 2 ' 



(71) 



Here the LM subscript in the boundary values CjL have A similar derivation provides the matrix elements Q r , 
been suppressed. Inserting Eq. (Ffll into Eq. ((70) yields given in Eq. (46). 



Eq. (44), 






P = 

J m,n. 



dfci? 2 



for n 7^ to (72) 



-m °n/ 



n ^V " ) d k R 2 ' 

k K 



